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Abstract 

The three-body problem in one-dimension with a repulsive inverse square potential 
between every pair was solved by Calogero. Here, the known results of supersymmetric 
quantum mechanics are used to propose a number of new three-body potentials which 
can be solved algebraically. Analytic expressions for the eigenspectrum and the eigen- 
functions are given with and without confinement. 
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I. Introduction 

Long back, in a classic paper (Calogero 1969), Calogero had given the complete solu- 
tion of the Schrodinger equation for three particles in one dimension, interacting pairwise 
by two-body harmonic and inverse-square potentials. Later, Wolfes (Wolfes 1974) had 
used Calogero's method to obtain analytical solutions of the same problem in the presence 
of an added three-body potential of a special form. Attention soon shifted to the exact 
solutions of the many-body problem (Sutherland 1971, Calogero 1971) and the general 
question of integrability. A list of solvable pair potentials for the many-body problem is 
given in the review by Olshanetsky and Perelomov (Olshanetsky and Perelomov 1981, 
1983), where other references will be found. Recently, there has been a renewed interest 
in the one-dimensional many-body problem of the Calogero and the Sutherland types 
(Polychronakos 1992), and their applications to the physics of spin chains (Haldane,1988, 
Sliastry 1988, Frahm 1993) .In this paper, we give some new potentials that have exact 
solutions of the three-body problem in one-dimension. Examples of confined spectra of 
the Calogero- Wolfes type, as well as other cases with both bound and scattering states 
present, are given. 

Three particles in one dimension, after the cetre-of-mass motion is eliminated, have 
two independent degrees of freedom. This may therefore be mapped on to a one-body 
problem in two- space dimensions as was done by Calogero. For the type of potentials 
used, the radial and angular parts separated, and the eigenvalue problem rested essen- 
tially on the availability of an analytical solution for the angular part. The same method 
is followed here, but advantage is also taken of supersymmetric quantum mechanics . The 
examples with solutions in a closed form that are presented here are found to be shape 
invariant supersymmetric potentials (Gendenshtein 1983, Dutt et al 1988, Khare et al 
1991) in the angular variable. Bound state algebraic solutions with harmonic confinement 
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are given in Section II, while scattering is discussed in Section III. The harmonic confine- 
ment of Section II may be replaced by a ^-type potential, and this is briefly discussed in 
the final Section. 



II. New examples with harmonic confinement 

Calogero had given the solution for the following potential for the three-body prob- 
lem : 

Vc = a;V8 J^ixi - Xjf + g ^(x^ - Xj)-^ (1) 

i<j i<j 

where g > —1/2 to avoid a collapse of the system. Wolfes (Wolfes 1974) showed that a 
three-body potential 

Vw = fiixi + X2- 2x3)~^ + ix2 + X3- 2xi)-^ + (x^ + xi - 2x2)~^] (2) 

is also solvable when it is added to Vc, with or without the pair- wise centrifugal term. The 
last two terms on the right-hand side of Eq.(2) are just cyclic permutations of the first. 
Henceforth, such terms occuring in other potentials would be referred to as 'cyclic terms 
'. In this section, we give more examples of three-body potentials that may be solved 
exactly. Our first solvable example is the three-body potential of the form 

which is added to Vc of Eq.(l). In the above equation, 

= ^[{xi - X2f + {X2 - x^f + {x^ - x^f]. (4) 
To see why potentials given by Eqs(l)-(3) are solvable, define the Jacobi coordinates 



R = ^{Xl +X2+ X3), 



_ (a^i - X2) _ {xi + X2- 2x3) 

After elimination of the centre-of-mass part from the Hamiltonian, only the x- and the 
y-degrees of freedom remain, which may be mapped onto the polar coordinates 

X = r sin 0, y = r cos 0, (6) 

where we follow the notation of Calogero. Obviously, the variables r, have ranges < 
r < 00, and < < 27r.It is straight- forward to show that 

(a^i — x-2) = ^f2 r sin 0, 
{x2 -X3) = V2r sin(0 + 27r/3), 

(xs-xi) = V2r sm{(t> + 47T/3). (7) 

The Schrodinger equation separates cleanly in the radial and angular variables, and the 
wave function may be written as 

■^ni{r,(j)) = Rni{r)Fi{(j)). (8) 

The radial wave function obeys the equation (7i = 1, 2M — 1). 

+ l^'r' + Bf/r')Rr^l = E^iR^i{r), (9) 
dr^ r dr 8 

where -B^^is the eigenvalue of the Schrodinger equation in the angular variable. The 
eigenvalues E^i are given by 

Eni = 7372 ^(2n + B1 + I). n = 0, 1, 2...; 1 = 0,1, 2... (10) 

The form of Rni is given by (Calogero 1969) 

Rni = {rf^ exp[^(y372)a;r2]Lf'[i(y372)a;r2], (11) 



whereB/ > O.To examine the angular part of the eigenfunction F/(0),take the po- 
tential Vc + Vi-Then, in the variable (j), the Schrodinger equation is 




(12) 

Calogero used the identity 
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cosec^[0 + 2(m - l)n/3] = 9cosec^(3(^). (13) 

m=l 

A similar identity may be used for cot (f), 

3 

cot[(p + 2(m - l)7r/3] = 3cot(30). (14) 

m=l 

Using Eqs(13) and (14), (12) reduces to 

[-^ + ^ ^cosec2(30) + ^ /i cot(30)] Fi = Bf Fi{^), (15) 

whose solution will be given presently. 

The solution of Eq(15) and some other examples to be given are found by algebraic 
methods of supersymmetric quantum mechanics (Infeld and Hull 1951, Gendenshtein 
1983). Define a "super potential" W{a'K)^ where a is a constant and the angular 
variable. Then the supersymmetric partners V- and Vj^ are {W'^ — W) and {W"^ + W) 
respectively,the prime denoting a differentiation with respect to 0.The potential V- has 
a zero energy ground state with the wave function Fq given by eyi^{— j'^ Wdcj)). In 
the literature (Dutt et al 1987, Levai 1989, Khare et al 1991 ), various superpotentials 
W are given for which V-,Vj^ are exactly solvable. In these examples, V_ and Kf are 
"shape invariant", meaning thereby that V+ is identical in functional form to V- with 
an appropriate change in one or more parameters. The potential in Eq(15) is obtained 
if we choose 

W = -AcotZ^-C/A, (^>0), (16) 
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where C,A are parameters independent of 0.Then 

= A{A ^ 3)cosec230 + 2Ccot 30 + {C'^/A^ - A^). (17) 
In this example, V± are shape invariant potentials, i.e. they are related as 

V+{A, C, 0) = V_{A + 3, C, 0) + (v4 + 3)2 - C'^/{A + 3f + C^/A^ - A^. (18) 
This in turn implies that the energy eigenvalues of V- are given by 

= {A + 3lf - Cy{A + 3lf -A^ + C^/A^. (19) 

From this, it immediately follows that for the potential 

9 o 9 

V = -g cosec^3(f) + - /icot 3(j) (20) 

occuring in Eq(15), the eigenvalue spectrum is 

Bf = 9(/ + a + 1/2)2 _ ^^2/^; ^ ^ ^ 1/2)2, ^21) 

where 

a = 1/2(1 + 2^)^/2 _ (-22) 

As in Calogero, g > —1/2 for meaningful solutions. The unnormalized ground state wave 
function is given by 

Fo(0) = exp[- Wd(t>\ = {sm3(t>r+-^ exp[^/i0/(a + 1/2)], (23) 

valid in the rangeO < 30 < tt. For distiguishable particles, a given value of defines 
a specific ordering. For < < 7r/3, Eq(7) implies xi > X2 > a;3,and other ranges of 
correspond to different orderings (Calogero 1969, Wolfes 1974). For singular repulsive 
potentials, crossing is not allowed, and Fo(0) of Eq.(23) is zero outside < < 7r/3. 
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Following Calagero, the wave function for the other ranges may be constructed. For indis- 
tiguishable particles, similarly, symmetrised or antisymmetrised wave functions may be 
constructed and this will not be repeated here. From the table given by Levai (Levai 
1989) the general solution of Eq.(15) is given by 



= exp(-z^O(sin30)'exp[^4^]p/ '^'^ ' (z cot 30). (24) 



In the above equation, P^'l^ is the Jacobi polynomial of the arguement(icot30). Another 
exactly solvable potential is given by 



i<j 

+ 3g [{xi + X2 — 2x3)"^ + cyclic terms ] 

Although this looks very different from the previous example of {Vc + Vi), it is simply 
obtained by the transformation^ ^ (j) + n/2. A little algebra shows that Eq.(25) reduces 
to 

Son q 9. , 2(m — l)7r, 3 fi v-^ , , 2(m — 

= 8^'^' + 2^ E + 3 ] - 2 ^ E + 3 ] 

m=l m=l 

= ^uj\^ + 9--^sec^30 - 9;^ tan 3(f) . (26) 
8 2r^ 2r^ 

The angular part of V2 is identical to Eq.(20) if ^ + 7r/2. Due to this shift, the 
allowed range of will divide into sectors — 7r/6 < </> < 7r/6, 7r/6 < </> < 7r/2, etc 
(Wolfes 1974). Otherwise, the solution Fi{(f)) is the same as Eq.(24) with (p + tv/2, 
with Bf stiU given by Eq.(21). There are more examples to be found by exploiting 
super symmetric quantum mechanics. A potential that is exactly solvable, but is different 
from the example given by Wolfes (Wolfes 1974) is 

Vs = \u^'Y.^x,-x,f + gY,{x,-x,)-'-^[ ^'')^''' -'^^'^ + cyclic terms]. (27) 
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Using relations (7) and (13), and the identity, 

3 

cos(0 + 2(m- l)7r/3)/sin^(0 + 2(m- l)7r/3) = 9cos 30/sin^30, (28) 

m=l 

V3 reduces to the form 

V3 = -a;2^2 ^ _ ^cosec^30 - - ^cot 30 cosec 30. (29) 
8 2 2 

As before, the Schrodinger equation is separable in r and 0, with the radial solution 
being the same as in Eq.(ll), and the angular part reduces to 

[-X^ + %cosec230- ^/3Cot30cosec3(/>]F,((/>) = Bf Fi{(P). (30) 
a(f) 2 2 

This form of V_(30) is obtained by taking the superpotential ( Dutt etal. 1988 )W to be 

= cot 30 + C cosec 30 (31), 

with < 30 < TT, ^ > 0, and C < ^.The resulting F_ is then given by 

V- = (A^ + C^- 3^)cosec^30 - (2A - 3)Ccot 30 cosec30 - A^, (32) 

which is the same as in Eq (30) apart from the constant^^ . The eigenvalues Bf are 
then independent of C, and given by 

Bf = {A + Zlf, Z = 0,1,2, . . (33) 

These are of the same form as obtained for the pair-wise inverse square potential. The 
unnormalised angular wave function F/(0) is given by 

FM) - (sin^)ncos^)^Pr^/^'"-^/^ (COS30) , (34) 

where a = (A + i?)/3, (3 = {A — B) J?,. Note that for the wave function Fi to vanish at 
30 = and at 30 = tt, a and (3 should both be > O.The range of the variable may be 
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extended following Calogero ( Calogero 1969 ). A different form of a solvable potential 
may again be obtained from Eq.(29) by transforming — > + 7r/2. Prom Eq.(6), this 
changes a; — > j/ and j/ — > — x.The resulting form of the solvable potential is given by 



V4 = -w^ ~ ^j)^ + 3^ [ {xi + X2 — X3) ^ + cyclic terms] 

H — ^ — [ 7 — — — ^l! + cyclic terms ] . (35) 
3a/2 r ^{xi + X2- 2x3)^ 

The reader may obtain the eigenvalues and the eigenfunctions of this potential (valid for 
the sectors — 7r/6 < (j) < 7r/6, 7r/6 < (j) < 7r/2 etc.) from the solutions of V3, and the 
details would not be given here. 

There are other types of three-body potentials that may be constructed whose an- 
gular part is solvable, but is not of the supersymmetric shape-invariant form as before. In 
such examples, however, although the eigenfunctions Fi{(j)) are expressible in terms of 
well-known functions, the eigenspectrum of Bf is not found to be in a simple closed 
form. For example,keeping in mind the identity, 

-{x\ - X2){x2 -X3){x3 - xi) = ^ sin 30, (36) 

v2 

a three-body potential of the form ^ sin^3(/) may be constructed. Combining this with 
the pair potential of Calogero that yielded l/(r^ sin^30), the angular part of the three- 
body Schrodinger equation may be written in the form 



—J- + [Bi - acos60 - — ^ 

sm 3(f) 



+ [Bf - acos60 - — — ] = 0. (37) 



It is not difficult to show that this may be reduced to the well-known diflFerential equation 
(Erdelyi 1955) for spheroidal wave functions (for6>a — 1/2), whose eigenvalues, even 
though known, are not expressible in a simple closed form. 
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Till now, we have been discussing only those cases where the Schrodinger equation 
in the variables r and (f) are separable, and exactly solvable. The three-body potential 
given by Eq.(36) is of special interest, however, even though its form does not allow 
separation of the coordinates r and 0. To be specific, consider the potential 

Vh = ^t^^ {Xi - Xj)'^ - X{Xi- X2)ix2 - X3){X3 - Xi). (38) 
i<j 

Transforming to polar coordinates using Eqs (7), this reduces to 

3 A 

Vh= -uj^r^ + sin 30, (39) 

which is the famous Henon-Heiles potential (Henon and Heiles 1964 ) for a particle in 
2-space dimensions. This potential is nonintegrable, and has been extensively studied in 
connection with chaos. Therefore Vh given by Eq(38) may be regarded as an example 
of a potential that gives rise to classical chaotic motion of three particles in one spatial 
dimension. Since it is well-known that classical periodic orbits have a close connection 
with the quantum density of states (Gutzwiller 1990), it will be of interest to map the 
two-dimensional periodic orbits of the Henon-Heiles potential (Brack et al. 1993) onto 
the one dimensional motion of the three particles. 

III. The three-body scattering problem 

The three-body scattering problem with the inverse square pair potential, 

i<j 

was solved by Marchioro (Marchioro 1969). When an additional three-body potential 
Vw (see Eq.(2)) is added to this, the problem is still solvable (Calogero and Marchioro 
1974). Both the classical and quantum solutions are strikingly simple in these problems. 
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Asymptotically, with the inverse square pair potential, the momenta of the scattered 
particles, are related to the incident ones pi by p^ = p^-i , {i = 1, 2, 3 ), and the total 
phase shift is independent of I. With the addition of the three-body potential Vw, the 
asymptotic momenta just change sign, and the phase shift still remains /-independent. We 
have also studied the scattering problem with the new potentials (cited in the last sec- 
tion ) after dropping the harmonic confinement. Even though the problems are still 
integrable, we find that the simplicity of the Calogero - Marchioro examples is lost. The 
scatterings with the new interactions cannot simply be expressed as an exchange between 
the incoming and asymptotic outgoing momenta. As an example, consider the classical 
scattering of three particles with the potential 

Vs=9j2i^i- ^j)-" - l^ll[i^L±^l^lpl + cyclic terms ], (40) 

which is the same as V3 given by Eq.(27), but without the harmonic confinement. Because 
the variables r, are separable, we may write, following Marchioro's notation ( Marchioro 
1969 ), 



1 o 52 
2 



pI + — = E, (41) 



1 9 o 9 

-pI + -9 cosec^ 30 /s cot 30 cosec 30 = 5^ . (42) 

Here E is the total energy and B the angular constant of motion.It is straight-forward 
to perform the integrations. Writing 

h 9-^3 ^ 18^ 

we obtain 

S2 - 



rit) = [2Eit-to)' + (44) 



and 



COS30 =1- lV4-c + b^sin [sin-^( ^/'"'^f° ) + 3t^n-\V2 ^ {t - to))] . (45) 
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In the above, = (f)o a,t t = to-Let (f) = (f)i for (t — to) — > oo, and (f) = (f)f for 
{t - to) oo.Then 

„, b 1 /- — r • ,b — 2cos3(bQ 

cos30i = - - -\/4-c + 6^cos[sin ( , ^ =r)] , 

^ ^ V 4 — c + 0^ 

and 

^1 r. r . 1 ,6 — 2 cos 3^0 NT /.^x 

COS30/ ^4 - c + 6^ cos[sin-^( ^ )] . (46) 

Therefore only for 6 = 0, i.e., for fs — ,does cos3(/)i = — cos3(/)/, leading to 3(/)i = 
(±30/ + tt). If one considers the sector for which < (f) < 7r/3, then 0j = —0/ + 7r/3 
for 6 = 0. It was this peculiarity that yielded the simple relationships between the initial 
and the final asymptotic momenta of the particles. But for /s ^ 0, no such relationship 
exists, and the simplicity in the scattering process is lost. For the same example, similar 
complications arise in the quantum treatment of the scattering problem with /s 7^ 0. In 
this case the angular wave function given by Eq.(34) does not possess a simple symmetry 
property for ^ — 7r/3 when (3. For other examples like Vi (given by Eq.(3)) plus 
the inverse square pair potential, the situation is even more complicated, and will not be 
discussed further 

IV. Discussion 

The shape-invariant potentials in the angular variable were combined with a har- 
monic confinement to obtain exact solutions in Section 2. This led to a discrete eigenvalue 
spectrum. The harmonic confinement may be replaced by an attractive i-type interac- 
tion,giving rise to both discrete and continuous energies. Exact solutions for all the 
shape-invariant potentials described earlier can again be obtained algebraically. As one 
example, take the potential 

V = . = + g y{Xi — Xj) , 

{Xi - X2)'^ + {X2 - Xs)'^ + {X3 - Xi)'^ 
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= -r + J -T-:^ ■ (47) 



a 9 q 

h — 

r 2 r^sin^S^ 

First consider the classical problem. Using the notation of Eqs. (41) and (42), we obtain, 



2B^ I 4EB^ A COS30 

=l + Jl + — cos-cos~^( , ) ■ (48) 

ar V a2 ^3 ^ ^1 - 9 g/{2 B^)'^ ^ ' 

For (7 = 0, and i? < 0, Eq.(48) reduces to the periodic orbit of an ellipse in the polar 
coordinates : 

- = 1 + e cos (/) , (49) 
r 



9 / 4l £/ 1 B'^ 

where p = ^^is half the latus rectum, and e = y 1 — is the eccentricity of the 

ellipse(Landau and Lifshitz 1960). For g = 0, crossing between any pair of particles is 
allowed, and the periodic orbit given by Eq.(49) may be mapped on to the motion of 
three particles along a line. On the other hand, for g ^ 0,Eq.(48) does not reduce to 
a closed orbit for the bound problem. For the special case when g ^ between one 
pair, and between the other two pairs, Eq.(48) still reduces to the closed orbit form of 
Eq.(49). However, crossing between the interacting pair is not allowed, and the mapping 
in each sector has to be done accordingly. The Schrodinger equation with the potential 
V of Eq.(47) is easy to solve. Writing the wave function i/^niir, (/)) — Rnio{r)Fi{(f)) as 
before, the radial part obeys the standard Coulomb-type equation with the bound-state 
wave function 

Rni{r) = r^' exp{-^/\E;^\r) '(2r/ro) , (50) 

where ro = 2(S/ + n + |)/q!, and L^^' is the Laguerre polynomial. The corresponding 
eigenvalues are 

^-'= , p , 1x2 ' ^ = 0,1,2,.. Z = 0,1,2... (51) 

4(n + Bi + 

As before, the angular constants Bis are the eigenvalues of the equation 

l-^m)= B^Fi{ct>). (52) 



2 sin^ 30 
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The unnormalised solutionsF/'s are expressed in terms of the Gegenbauer polynomials 
(Calogero 1969) 

Fi = (sin30)"+^ C"+^ (cos30) , < < ^ . (53) 

Here the constant a is given by Eq.(22), and the wave function vanishes outside the range 
xi > X2 > X3. The eigenvalues Bf are the same as in Eq.(21) with /i = 0, so 

Bi = 3il + a+^). (54) 

Similarly, for the continuum states, the scattering problem for the radial Coulomb part 
may be solved in the standard fashion (Landau and Lifshitz 1960), with the phase shift 
Si given by the usual Coulomb expression,but with I replaced by Bi. 

The three-body problem in higher dimensions may also be mapped onto a one- 
body problem. For example, the three-body problem in 2-space dimensions has been 
reduced to the 4-dimensional hyperspherical coordinates (Kilpatrick and Larsen 1987) 
of one particle. This method has been used to obtain the spectrum of three-anyons in 
an oscillator potential (Khare and McCabe 1991, Law et al. 1992). The technique of 
supersymmetric quantum mechanics, so aptly adapted for the one-dimensional three- 
body problem, has not been generalised for higher dimensions , and hence algebraic 
solutions cannot be obtained in such cases. 
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